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Abstract 



We present a comprehensive phenomenological analysis of a vast amount of data from neutrino 
^r) ' flavor oscillation and non-oscillation searches, performed within the standard scenario with three 

massive and mixed neutrinos, and with particular attention to subleading effects. The detailed 
results discussed in this review represent a state-of-the-art, accurate and up-to-date (as of August 
2005) estimate of the three-neutrino mass-mixing parameters. 
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^ ! 1 Introduction 

in 

O | Neutrinos provide, on a macroscopic scale, the realization of two key concepts of quantum mechanics: 
^ linear superposition of states and noncommuting operators. In fact, there is compelling experimental 
evidence [Tj that the three known neutrino states with definite flavor (V e , z/ M and v T ) are linear combi- 
Ph. nations of states with definite mass z/, (i = 1,2,3), and that the Hamiltonian of neutrino propagation 
in vacuum [2] and matter [3J HJ E] does not commute with flavor. The effects of flavor nonconservation 
j>! | ( "oscillations" ) take place on macroscopic distances, for typical ultrarelativistic neutrinos. The evidence 
for such effects comes from a series of experiments performed during about four decades of research 
with very different neutrino beams and detection techniques: the solar neutrino [0] experiments Home- 
! stake [T], Kamiokande g], SAGE 0, GALLEX-GNO [10111111121, Super-Kamiokande (SK) [HUH] and 
Sudbury Neutrino Observatory (SNO) 011011121; the long-baseline reactor neutrino [TH] experiment 
KamLAND [T§ll2"Ul l2"T]: the atmospheric neutrino [22] experiments Kamiokande [20], Super-Kamiokande 
[2H ESI EE] , MACRO (2Z|, and Soudan-2 j2E]; and the long-baseline accelerator neutrino [2H| experiment 
KEK-to-Kamioka (K2K) [SHIISI]. 

Together with the null results from the CHOOZ [32] (and Palo Verde [00]) short-baseline reactor 
experiments, the above oscillation data provide stringent constraints on the basic parameters governing 
the quantum aspects of neutrino propagation, namely, the superposition coefficients between flavor and 
mass states (i.e., the neutrino mixing matrix), the energy levels of the Hamiltonian in vacuum (i.e., 
the splittings between squared neutrino masses), and the analogous levels in matter (i.e., the neutrino 
interaction energies). The energy levels at rest (i.e., the absolute neutrino masses) are being probed by 
different, non-oscillation searches: beta decay experiments [031 ESI EE] , neutrinoless double beta decay 
searches (GV2/3) (371 EE1 EH] , and precision cosmology [101 EU H21 1321 El] ■ Current non-oscillation data 
provide only upper limits on neutrino masses, except for the claim by part of the Heidelberg-Moscow 
experimental collaboration ^Sl EOl IS], whose possible 0V2/3 signal would imply a lower bound on 
neutrino masses. 
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A highly nontrivial result emerging from these different neutrino data sample is their consistency, 
at a very detailed level, with the simplest extension of the standard electroweak model needed to ac- 
commodate nonzero neutrino masses and mixings, namely, with a scenario where the three known 
flavor states ^ e , ftT are mixed with only three mass states ^1,2,3, no other states or new neutrino in- 
teractions being needed. This "standard three- neutrino framework" (as recently reviewed, e.g., in 
[HI EHl EIIJ EH E21 EH1 EH ESI EH I37j ) appears thus as a new paradigm of particle and astroparticle 
physics, which will be tested, refined, and possibly challenged by a series of new, more sensitive exper- 
iments planned for the next few years or even for the next decades IHH]. The first challenge might 
actually come very soon from the running MiniBooNE experiment jS0j , which is probing the only piece 
of data at variance with the standard three-neutrino framework, namely, the controversial result of the 
Liquid Scintillator Neutrino Experiment (LSND) [p"!] . 

In this review we focus on the current status of the standard three-neutrino framework and on the 
neutrino mass and mixing parameters which characterize it, as derived from a comprehensive, state-of- 
the-art analysis of a large amount of oscillation and nonoscillation neutrino data (as available in August 
2005). All the results and figures shown in the review are either new or updated or improved in various 
ways, with respect to our previous publications in the field of neutrino phenomenology. In this sense 
we have tried to be as complete as possible, so as to present a self-consistent overview of the current 
status of the three-neutrino mass-mixing parameters. While we aimed at obtaining technically accurate 
and complete results, we have not aimed at being bibliographically complete; we refer the reader to 

HU 1221 IS EOl E21 EH E3 EH1 EH EH1 COl IZH H21 113 EH 1^1 113 HZl IS EH for an incomplete 

list of excellent reviews with rich bibliographies of old and recent neutrino papers. 



2 Notation 



While for quark mixing a standard notation and parametrization has emerged PQ, this is not (yet) 
the case for (some) neutrino mass and mixing parameters. In this section we define and motivate the 
conventions used hereafter. 



2.1 Mixing angles and CP-violating phase 

At the lagrangian level, the left-handed neutrino fields with definite flavors v a L (a = e, /i, r) are assumed 
to be linear superpositions of the neutrino fields with definite masses vn (i = 1, 2, 3), through a unitary 
complex matrix U a f. 

3 

V<xL = J2U ai V iL . (1) 

1=1 

This convention implies [H21IHD1 that one-particle neutrino states \v) are instead related by U* (see, e.g., 

fflEU), 

Wa) = J2KM)- (2) 

i=l 

A common parameterization [HI] for the matrix U is: 

U = 23 T 5 13 Tl0 12 , (3) 

where the Oy's are real Euler rotations with angles e [0,7r/2] |B2], while T$ embeds a CP-violating 
phase 5 G [0, 2%], 

r 5 = diag(l,l,e +l5 ) . (4) 

Notice that the above definitions imply det(U) = +1, which may be a useful property in some theoretical 
contexts 
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By considering TsOi^Tl as a single (complex) rotation, this parametrization coincides with the one 
recommended (together with Eq. (J2J)) in the Review of Particle Properties PP, 



U 



1 
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V 01 



(5) 



where Cy = cos 6^- and s y - = sin%. 

Other conventions sometimes used in the literature involve U instead of U* in Eq. (J2J), or —5 instead 
of +5 in Eq. or only one CP- violating factor (either T$ or T\, not both) in Eq. @, or a combination 
of the above. In our opinion such alternatives, although legitimate, do not bring particular advantages 
over the above convention. 

For the sake of simplicity, the phase 5 will not be considered in full generality in this work. Numerical 
examples will refer, when needed, only to the two inequivalent CP-conserving cases, namely, e lS = ±1. 
In these two cases, the mixing matrix takes a real form Uqp, 



C13C12 

tA:p = I -S12C23 =F S23S13C12 

<S23 S 12 T S13C23C12 



S12C13 ±Si3 
C23C12 =F S23S13S12 S23C13 1 • 
— S23C12 =F S13S12C23 C23C 13 



(6) 



where the upper (lower) sign refers to 5 = (5 = n). The two cases are formally related by the 
replacement S13 — > — S13. In any case, CP violation effects do not affect at all solar and reactor 
oscillation searches, where the indistinguishability of and v T in the final state allows to rotate away 
both the angle #23 and the CP phase S from the parameter space, even in the presence of matter effects 
(see, e.g., jHlj). 



2.2 Masses, splittings and hierarchies 

The current neutrino phenomenology implies that the three-neutrino mass spectrum {w^}i=i,2,3 is 
formed by a "doublet" of relatively close states and by a third "lone" neutrino state, which may be 
either heavier than the doublet ("normal hierarchy," NH) or lighter ("inverted hierarchy," IH). 1 In the 
most frequently adopted labeling of such states, the lightest (heaviest) neutrino in the doublet is called 
u i (^2), so that their squared mass difference is 

5m 2 = m\ — m\ > (7) 

by convention. The lone state is then labeled as u 3 , and the physical sign of m\ — m\ 2 distinguishes 
NH from IH. 2 

Very often, the second independent squared mass difference Am 2 is taken to be either m\ — m\ or 

"3 



m\ — m\. However, these two definitions may not be completely satisfactory in both hierarchies. In 



fact, in passing from NH to IH, the difference m\ — ml not only changes its sign, but also changes 
from being the largest squared mass gap to being the next-to-largest gap (while the opposite happens 
for m 2 — m 2 ,). Whenever terms of 0(5m 2 / Am 2 ) are relevant, this fact makes somewhat tricky the 
comparison of results obtained in different hierarchies. For such reason, we prefer to define Am 2 as 



Am 2 



m 3 - 



2 

1 In this context, "hierarchy" does not refers to neutrino masses but only to mass differences. In particular, it is not 
excluded that such differences can be much smaller than the masses themselves — a scenario often indicated as "degenerate 
mass spectrum." 

2 Another convention, sometimes used in the literature, labels the states so that mi < iri2 < in both NH and IH. 
In this case, however, the mixing angles 6ij have a different meaning in NH and IH. 
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so that the two hierarchies are simply related by the transformation +Am 2 
next-to-largest squared mass gaps are given Am 2 ± 5m 2 / 2 in both cases. 
More precisely, the squared mass matrix 

M 2 = diag(m 2 , m 2 ,, m 2 ) 



—Am 2 . The largest and 



(9) 



reads, in our conventions, 



M 



mf + m 2 . 



( Sin 2 Stti 2 
1 + diag , H , ±Am 2 



(10) 



where the upper (lower) sign refers to normal (inverted) hierarchy. 

In the previous equation, the term proportional to the unit matrix 1 is irrelevant in neutrino os- 
cillations, while it matters in observables sensitive to the absolute neutrino mass scale, such as in 
/3-decay and precision cosmology. In particular, we remind that /3-decay experiments are sensitive to 
the so-called effective electron neutrino mass nip, 



E 



|2 2 



2 2 *l \ jl- — _l . _ 
C 13 C 12 m l + C 13 S 12 m 2 + S 13 m . 



2 „2 



:ir 



as far as the single V{ mass states are not experimentally resolvable [8*7| . On the other hand, precision 
cosmology is sensitive, to a good approximation (up to small hierarchy-dependent effects which may 
become important in next-generation precision measurements 88 ) to the sum of neutrino masses £ 



£ = nil + m 2 + m 3 ■ 



(12) 



2.3 Majorana phases 

If neutrinos are indistinguishable from their antiparticles (i.e., if they are Majorana rather than Dirac 
neutrinos), the mixing matrix U acquires a (diagonal) extra factor [0211111113] 



M 



(13) 



which is parametrized in various ways in the literature. In particular, within the Review of Particle 
Properties, two different conventions are used |53*1 192j. We adopt the one in [02], which — after a slight 
change in notation — reads: 

U M = diag (l, e^ 2 , e5 {<fe+2<5) ) , (14) 

4>2 and 03 being unknown Majorana phases. The "advantage" of this convention is that, in the expression 
of the effective Majorana mass mpp probed in neutrinoless double beta decay (Oz/2/?) experiments 
[231021, the CP- violating phase 5 is formally absent: 



c i3 c i2 m i + Ci 3 s 12 m 2 e l<?2 + s l3 m 3 e l 



(15) 



2.4 Matter effects 

In the flavor basis, the hamiltonian of ultrarelativistic (m, p) neutrino propagation in matter reads 



H 



IE 



UM 2 U^ + K 



MSW 



(16) 
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up to an irrelevant momentum term pi which, acting as a zero-point energy, produces only an unob- 
servable overall phase in flavor oscillation phenomena. In the above equation, Vmsw = diag( V, 0, 0) 
is the Mikheyev-Smirnov-Wolfenstein (MSW) term [3| embedding the interaction energy difference (or 
"neutrino potential"), 

V(x) = V2G F N e (x) , (17) 

E being the neutrino energy, and N e the electron density at the position x. For antineutrinos, one has 
to replace U — > U* and V — *■ — V. We shall also use an auxiliary variable with the dimensions of a 
squared mass [55] , 

A(x) = 2EV = 2V2G F N e (x)E . (18) 

Matter effects are definitely important when one squared mass difference (either 5m 2 or Am 2 ) is of the 
same order of magnitude as A(x). 

When needed, the eigenvalues of H in matter will be denoted as m 2 /2E, and the diagonalizing 
matrix as U (with rotation angles 9^): 

H = 7^UM 2 W . (19) 

The eigenvalue labeling is fixed by the condition mf — > m\ for A[x) — > 0. The parameters rrtj and 8ij 
are often called "effective" neutrino masses and mixing angles in matter. 

We remind that, in the absence of matter effect, and within the two CP-conserving cases (e lS = 
±1 — ► U — U*), the (vacuum) flavor oscillation probability P a p = P(u a — ► vp) takes the form [HB] 

(jji^ — Tn 2 \ 
\ E J Lj , (20) 

where L is the neutrino pathlength. The same functional form is retained in matter with constant 
density, but with mass- mixing parameters (6ij,m 2 — m 2 ) replaced by their effective values in matter 



lij,m 2 -m|) 



' ~ 2 ~ 2 

mf — m- 



P^T = <W - 4 E U ai U aj U pi U Pj sin 2 ^ ^ E J L ] . (21) 



For non-constant matter density, P a/ 3 cannot be generally cast in compact form and may require nu- 
merical evaluation, although a number of analytical approximations can be found in the literature for 
specific classes of density profiles. 



2.5 Conventions on confidence level contours 

In this work, the constraints on the neutrino oscillation parameters have been obtained by fitting 
accurate theoretical predictions to a large set of experimental data, through either least-square or 
maximum-likelihood methods. In both cases, parameter estimations reduce to finding the minimum of 
a x 2 function (see the Appendix) and to tracing iso-Ax 2 contours around it. 

Hereafter, we adopt the convention used in pQ and call "region allowed at ncr" the subset of the 
parameter space obeying the inequality 

A X 2 < n 2 . (22) 

The projection of such allowed region onto each single parameter provides the na bound on such 
parameter. In particular, we shall also directly use the relation V 'Ax 2 = n to derive allowed parameter 
ranges at n standard deviations. 
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3 Solar neutrinos and KamLAND 



In this section we present an updated analysis of the constraints on the mass-mixing parameters placed 
by oscillation searches with solar neutrino detectors and long-baseline reactors (KamLAND) in the 
parameter space (6m 2 , sin 2 # 12 , sin 2 13 ). We start with the limiting case 6> 13 — > 0, and discuss in detail 
the bounds on (5m 2 , #i 2 ). We also discuss the current evidence for the occurrence of matter effects in the 
Sun, and then describe some details of the statistical analysis. We conclude the section by discussing 
the more general case with #13 unconstrained. 

Some technical remarks are in order. The latest KamLAND results [201 EI] are analyzed through 
a maximum-likelihood approach including the event-by-event energy spectrum [HH]. Here we do not 
include the additional time information available in [2lj which, as discussed in [53], does not improve 
significantly the bounds on the oscillation parameters. Solar neutrino data are analyzed through the 
pull method discussed in [94]. With respect to [§3|, Chlorine j?j and Super-Kamiokande data [T3] 
are unchanged, while Gallium results have been updated [3 HH E2]- 111 addition to the SNO-I (no 
salt) results jTHj already discussed in [03], we include in this work the complete SNO-II data (with 
salt) [T7j, namely, day and night charged- current (CC) spectra (17+17 bins), and global day and night 
neutral-current (NC) and elastic-scattering (ES) event rates (2+2 bins), together with 16 new sources of 
correlated systematic errors affecting the theoretical predictions 17J. Correlations of statistical errors 
(treated as in [HE]) in SNO-II data [H] are also included. Some of the SNO systematics are highly 
asymmetrical and even one-sided [17], and their statistical treatment is not obvious. We have chosen 
to apply the prescription proposed in jHE] to deal with combinations of asymmetric errors: for each i-th 
pair of asymmetric errors (af,a~) affecting a theoretical quantity R, we apply the pull method 94 
to the shifted theoretical quantity R + ARi with symmetric errors ±<jj, where 2Ai?j = af — o~ and 
2(jj = af + er^ [96 . Care must be taken to account for relative bin-to-bin error signs. We understand 
that the SNO approach to asymmetric errors (not explicitly described in |17j ) is different from ours 
[97]; this fact might account for some differences in our allowed regions, which appear to be somewhat 
more conservative at high-5m 2 values, as compared with those in J7| . Finally, the input standard solar 
model (SSM) used in this work is the one developed by Bahcall and Serenelli (BS) in [§B1 HUH HOUj by 
using a new input (Opacity Project, OP) for the opacity tables and older heavy-element abundances 
consistent with helioseismology [MI]. In this SSM (denoted as "BS05 (OP)" in [100]), the "metallicity" 
systematics, previously lumped into a single uncertainty, are now split into 9 element components. In 
total, our solar neutrino data analysis accounts for 119 observables [1 Chlorine + 2 Gallium (total 
rate and winter-summer asymmetry) + 44 SK + 34 SNO-I + 38 SNO-II] and 55 (partly correlated) 3 
systematic error sources. Further technical details can be found in the Appendix. 

3.1 Solar and KamLAND constraints (#13 = 0) 

For #13 = 0, electron neutrinos are a mixture of v\ and z/ 2 only. So, the parameter space relevant for 
solar Ue's and KamLAND u e 's reduces to the two variables governing the (fi,W) oscillations, namely, 
Sm 2 and # 12 . Trigonometric functions useful to plot #12 are either tan 2 # 12 in logarithmic scale or sin 2 # 12 
in linear scale; these choices graphically preserve octant symmetry (# 12 — > vr/2 — # 12 ) when applicable 
(e.g., in the limit of vacuum oscillations). 

Figure 1 shows the current solar neutrino constraints from separate data sets (Chlorine, Gallium, 
SK, SNO) at the 2a level, using the BS05 (OP) SSM input [100] . In each panel, we also superpose the 
small region allowed at 2cr around 5m 2 ~ few x 10~ 5 eV 2 and tan 2 # 12 ~ few x 10 _1 , which provides the 
solution to the solar neutrino problem [T2| at large mixing angle (LMA). In Fig. 1 one can appreciate 
that the global LMA solution completely overlaps with each of the regions separately allowed by the 

3 All the SSM sources of uncertainties are independent, with the exception of some of those concerning the "new" SSM 
metallicities, whose correlations are taken as recommended in |1(JU) . 
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different experimental data (at 2a), i.e., there is a strong consistency between different observations. 
The shape of the global solar LMA solution appears to be dominated by SNO and (to a lesser extent) by 
the SK experiment. Since both SK and SNO are sensitive to the high-energy tail of the solar neutrino 
spectrum (i.e., to the 8 B neutrino flux |lUUj ). and since the SNO determination of the 8 B flux is already 
a factor of two more accurate than the corresponding prediction in the BS05 (OP) standard solar model 
(see next Sec. 3.2), the shape of the global LMA solution in Fig. 1 is rather robust with respect to 
possible variations in the standard solar model input (including those related to the recent chemical 
controversy about the solar photospheric metallicity [0H1 EM D-OOJ). 

Notice that current solar neutrino data, by themselves, identify a unique (LMA) solution in Fig. 1; 
this was not the case only a few years ago (see, e.g., |10HI102| ). when at least another region at low Sm 2 
("LOW" solution) was allowed . From a test of hypothesis, we get that the current probability of 
the LOW solution is only Plow = 1-2 x 10 -3 . Former solutions in the vacuum oscillation regime (VAC) 
or at small mixing angle (SMA) (with acronyms taken, e.g., from |1U1[ 1102 j) are now characterized by 
exceedingly low probabilities (Pvac = 4.8 x 10~ 6 and Psma = 4.0 x 10~ 8 from current solar neutrino 
data). 

The LMA solution is heavily affected by solar matter (MSW) effects (see, e.g., |103U104j for a recent 
review of the LMA-MSW properties). Figure 2 shows the neutrino potential V(x) as a function of the 
normalized Sun radius x = R/R & , together with typical solar v e production regions (in arbitrary vertical 
scale), as taken from the BS05 (OP) model [23 COO.- F rom this figure one can easily derive that, for 
Sm 2 values in the LMA region, matter effects are definitely important (5m 2 ~ A(x)) for neutrinos with 
E ~ few MeV. More precisely, the LMA solar u e survival probability at the Earth [P ee = P[y e — > v e )\ 
reads ji QUI] 

Pee = 2 + \ cos20 12 O) cos2# 12 , (23) 

where 

cos 26 12 = cos26 12 -A(x)/5m 2 

(cos 29 12 - A(x)/5m 2 ) 2 + sin 2 26 12 



with cos 2#i2 slowly changing from its vacuum value (cos 29i 2 ) to its matter-dominated values (close to 
— 1) as E increases from sub-MeV to multi-MeV values. 

Figure 3 shows the energy profile of P ee , averaged over the production regions relevant to pp, 7 Be, 
and 8 B solar neutrinos, for representative LMA oscillation parameters. Also shown are the energy 
profiles of corresponding solar v e fluxes (in arbitrary vertical scale). The value of P ee decreases from 
its vacuum value (1 — 0.5 sin 2 26* 12 ) to its matter-dominated value (~ sin 2 # 12 ) as the energy increases. 
The vacuum-matter transition is faster for neutrinos produced in the inner regions of the Sun. In 
Fig. 3 we also show the small difference between day (D) and night (N) curves, due to matter effects 
in the Earth 4 (calculated, for definiteness, at the SNO latitude). The vacuum-matter transition is 
slightly slower during the night, due to the Earth regeneration effect (see |10fi| and references therein). 
Within current energy thresholds and experimental uncertainties, the vacuum-matter transition and 
the Earth regeneration effects have not been yet observed in the SK and SNO time-energy 
spectra. Nevertheless, as we shall see later, matter effects in the Sun must definitely occur to explain 
the data. 

Let us consider now the impact of KamLAND data. For typical LMA parameters, reactor v e are 
expected to have a relatively large oscillation amplitude (sin 2 2#i 2 ), as well as a sizable oscillation 
phase [5m 2 L/4:E ~ 0(1)] over long baselines (L ~ O(10 2 ) km). The v e disappearance signal observed 
in KamLAND 1201 nas n °t only confirmed the solar LMA solution but has greatly reduced its 
Sm 2 range [213 HZj, by observing a strong distortion in the energy spectrum (201 • Figure 4 shows 
the mass-mixing parameter regions separately allowed by the KamLAND total rate, by the energy 



4 The treatment of Earth matter effects in the present work is the same as in |1U5] but with eight density shells |94j. 
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spectrum shape, and by their combination, at the 1, 2, and 3a level, as obtained by our unbinned 
maximum-likelihood analysis [93J of the latest energy spectrum data j2Ij. The overall reactor neutrino 
disappearance (rate information) and its energy distribution (shape information) are highly consistent, 
the latter being dominant in the combination. At the 2a level, both the shape-only and the rate+shape 
analyses identify a single solution at 5m 2 ~ 8 x 10~ 5 eV 2 and large mixing; only at the 3a level two 
disconnected solutions appear at higher and lower values of 5m 2 . Notice the linear scale on both axes, 
and the reduction of the parameter space, as compared with Fig. 1. 

Figure 5 shows the regions separately allowed by all solar neutrino data and by KamLAND, both 
separately and in combination, at the 1, 2, and 3a level. The current solar LMA solution, as compared 
with results prior to complete SNO-II data (see, e.g., [TBI 1201 H07| . is slightly shifted toward larger 
values of sin 2 # 12 and allows higher values of 5m 2 . 5 This trend is substantially due to the larger value 
of the CC/NC ratio measured in the complete SNO II phase (0.34 |17j ) with respect to the previous 
central value (0.31 ^H])- We also find that the SNO-II CC spectral data ^7j contribute to allow 
slightly higher values of 5m 2 with respect to older results. The consistency of solar and reactor allowed 
regions is impressive, with a large overlap even at the la level, and with very close best-fit points. The 
solar+KamLAND combination eliminates the extra (KamLAND-only) solutions at high and low 5m 2 , 
and identifies a single allowed region characterized by the following 2a ranges: 

5m 2 = 7.92 x 10~ 5 eV 2 (1 ± 0.09) at ± 2a , (25) 

sin 2 12 = 0.314 (l±8:ii) at ± 2a . (26) 

The determination of these two parameters at O(10%) level represents one of the most remarkable 
successes of the last few years in neutrino physics. 

The 5m 2 uncertainty is currently dominated by the KamLAND observation of half-period of os- 
cillations (20] and can be improved with higher statistics [108.. The sin 2 #12 uncertainty is instead 
dominated by the SNO ratio of CC to NC events, which is a direct measurement of P ee at high energy: 
Rcc/R^c — Pee — sin 2 6' 12 . Figure 6 shows isolines of this ratio in the mass-mixing parameter space, 
which can be used as a guidance [109] IllOj to understand the effect of prospective SNO measurements 
on the sin 2 612 range. In the same figure we show isolines of the day-night asymmetry (Adn) of CC 
events in SNO, whose measurement could, in principle, help to reduce the 5m 2 uncertainty |109l IllOj: 
however, it is unlikely that the SNO errors can be reduced enough (< 1%) to clearly observe a day-night 
effect (see, e.g., |lllj ). 

3.2 Evidence for matter effects in the Sun 

As shown in Fig. 3, solar matter effects make P ee decrease from its vacuum value (> 1/2) to a matter- 
dominated value (< 1/2), for typical LMA parameters around the current best fit. Model-independent 
tests of the presence of matter effects are derived in the following, first by showing that SK and SNO 
data consistently indicate that P ee < 1/2 in their energy range, and secondly by showing that all 
solar -l-reactor data consistently indicate that the neutrino potential V(x) must be nonzero. 

As discussed in |112[ I113j . the normalized energy spectra of neutrinos which do interact in SK and 
in SNO (i.e., the SK and SNO "response functions" to the incoming 8 B neutrinos) can be equalized, 
to a good approximation, by choosing a proper SK energy threshold, for any given SNO threshold. 
The current best "equalization" of SK and SNO response functions is shown in Fig. 7. In this case, 
both SK and SNO are sensitive to the same energy-averaged survival probability, (P ee ). Moreover, the 
CC and NC event rates in SNO, together with the ES event rate in SK, overconstrain (P ee ) and the 
unoscillated 8 B solar neutrino flux $g in a completely model-independent way 6 (i.e., independently of 

5 Our current best-fit point for solar data only is at 8m? = 6.3 x 10 -5 eV 2 and sin 2 #12 = 0.314. 
6 For purely active (no sterile) neutrino flavor transitions. 
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the mass-mixing parameters and of the standard solar model) through the equations |112t I113| 

®B[(Pee)+ra(l-(P e e))} , (27) 
^B(Pee) , (28) 
®B , (29) 

where r a ~ 0.154 is the ratio of the energy-averaged ES cross-sections of and v e in SK. 

Figure 8 shows the current bounds at 2a on and (P ee ), as obtained by using the latest SNO CC 
and NC jT2j and SK ES [131 IHj event rates, both separately (bands) and in combination (2a elliptical 
regions). The dotted ellipse represent the combination of SNO NC and CC data; the addition of SK 
ES data — which are consistent with SNO NC and CC data — slightly increases the preferred value of 
$s (solid ellipse). In particular, the SNO+SK combination (dominated by SNO) provides the following 
ranges: 

$ B = 5.2t° i x lO^m-V 1 (±2ct) , (30) 

(Pee) = 0.34±SSS (±2(7) • (31) 

The above SNO+SK range for is consistent with the ±2(7 prediction of the BS05(OP) standard 
solar model [TUU] . $| SM = 5.7(1 ± 0.32) x 10 6 cm _2 s _1 , the difference in the central values (~ 10%) 
being not statistically significant, as also evident in Fig. 8. Notice that the SK+SNO data determine 
$s with an error a factor of 2 smaller than the SSM prediction. At the same time, the SK+SNO data 
constrain (P ee ) to be definitely less than 1/2 |114H115j . and in particular close to ~ 1/3, as predicted for 
high-energy 8 B neutrinos and LMA parameters (see Fig. 3). The model- independent SNO+SK analysis 
is thus fully consistent with the LMA-MSW expectations; removal of the MSW effect in the LMA region 
would give a prediction (P ee ) = 1 — 0.5sin 2 26> 12 > 1/2, inconsistently with the results in Fig. 8. 7 

One can perform, however, a more powerful test of the presence of the neutrino potential V(x), by 
artificially altering its magnitude through a free parameter Omsw (US IH5| lllfcij . 

V(x) -> ausw V(x) , (32) 

both in the Sun (relevant for solar neutrino oscillations) and in the Earth (relevant for both solar and 
reactor neutrino oscillations), and by renalyzing solar and KamLAND data with czmsw free. Testing 
matter effects amounts then to reject the case <2msw = (no effect) and to prove that cxmsw = 1 
(standard effect) is favored. In the analysis, we add CHOOZ reactor data, which help to exclude the 
appearance of spurious high-£m 2 solutions for cxmsw ^> 1 |114[ 11151 I116j . Figure 9 shows the results of 
a fit to all the current solar and reactor data in the parameter space (5m 2 , sin 2 # 12 , omsw); marginalized 
with respect to the first two parameters, in terms of the function (A^ 2 ) 1 / 2 = na. The preference for 
standard matter effects (gimsw = 1) is really impressive, and is currently even more pronounced then 
with previous data |116j . The hypothetical case of no matter effects (omsw = 0) is rejected at > 5a. 
Since V oc Gp, the results in Fig. 9 can not only be seen as a confirmation of matter effects, but can also 
be interpreted as an alternative "measurement" of the Fermi constant Gp through neutrino oscillations 
in matter, within a factor of ~ 2 uncertainty at 2a. 

3.3 Statistical checks 

We have seen that, globally, solar neutrino experiments agree with each other and with the KamLAND 
observation of reactor neutrino disappearance, that solar+KamLAND data identify a restricted range 
of LMA mass-mixing parameters, and that there is solid evidence for the associated matter effects in 

7 For (9 13 > 0, the no-MSW prediction would be slightly modified as (P ee ) > 0.5cf 3 + sf 3 = 0.46 (using the 3er upper 
limit sf 3 < 0.047 discussed later), still inconsistently with Fig. 8. 
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such range. However, it makes sense to look at the statistical consistency of the LMA best-fit solution 
in more detail, for at least two reasons: (1) the analysis involves a large number of observables and of 
systematics, some of which might deviate from the predictions without really altering the global fit; (2) 
the preferred shifts of some quantities might reveal something interesting. 

We remind that the solar neutrino analysis is performed through the so-called pull approach [94J, 
namely, by allowing shifts of each ra-th theoretical prediction i? n through independent systematic un- 
certainties C n k, 

Rn —* Rn = Rn + C nk i (33) 

k 

whose amplitudes are constrained through a quadratic penalty term. The shifted predictions -R n 's 
are then compared to the experimental values R e * p via the uncorrelated (mainly statistical) error com- 
ponents. The method can be generalized to include correlation of statistical jHS] and systematic |117j 
errors. The global x 2 function is then given by two terms, \ 2 — xlbs + X 2 ys y embedding the quadratic 
pulls of the observables (i.e., the deviations of theory vs experiment) and of the systematics (i.e., their 
offset with respect to zero). This method allow a detailed check of possible pathological deviations 
(pulls) of some quantities. (See also the Appendix.) 

Figure 10 shows the pulls of the 119 observables at the global (solar+KamLAND) best-fit point. 
From top to bottom, the observables include the Chlorine rate, the Gallium rate and its winter-summer 
asymmetry \11H\ IT2*] . the SK distribution in energy and zenith angle (44 bins) fSl Ej- the SNO-I 
(no-salt) CC spectrum in 17+17 day-night bins ^Sj, the SNO-II (with salt added) CC spectrum in 
17+17 day-night bins and the day-night values of the NC and ES rate The SNO-II data and 

their correlations are treated as prescribed in ^7j; for all the other observables we refer the reader to 
|94j . None of the pulls in Fig. 10 exceeds 3a, and their distribution, which is roughly gaussian, reveals 
nothing pathological. We conclude that none of the solar neutrino observables shows an anomalous or 
suspect deviation from the LMA best-fit predictions. 

Figure 11 shows the pulls of the 55 systematic errors which enter in the analysis. From top to bottom, 
they include 11 "old" standard solar systematics as in [01], 9 "new" SSM metallicity systematics [98J, 
the 8 B spectrum shape uncertainty |119j . 11 SK and 7 SNO-I systematics as in [H4] . and 16 "new" 
SNO-II systematics ^7]. All the offsets are small (< la), indicating that the allowance to shift the 
theoretical predictions R n through systematic uncertainties is only moderately exploited in the fit; in 
other words, there is no need to stretch the systematics beyond their stated la range to achieve a good 
fit. 

Finally, Fig. 12 shows a by-product of the pull approach, namely, the preferred shifts of the solar 
neutrino fluxes with respect to their central SSM values |94j . The ~ 10% downward shift of is 
consistent with the results in Fig. 8. The global fit also prefers a ~ 10% reduction of beryllium (Be) 
and CNO solar neutrino fluxes with respect to the BS05 (OP) prediction — an indication which may be 
of interest for future experiments directly sensitive to such fluxes |1201 1121j . Such preferred reductions 
are well within SSM uncertainties 100 . 

In conclusion, the detailed analysis of the LMA best-fit solution reveals a very good agreement 
between all single pieces of experimental and theoretical information in the solar neutrino analysis. No 
statistically alarming deviation is found. 

Concerning the statistical analysis of KamLAND data, the adopted maximum-likelihood approach 
93] involves only three systematic uncertainties, namely, two free background normalization factors a' 
and a" plus one constrained pull a for the energy scale offset (the overall rate normalization error being 
incorporated in the likelihood rate factor, see the Appendix and [22 )■ Therefore, our pull analysis for 
KamLAND involves a single parameter (a), which we find to be very small (a ~ 0.15o"q,) at the LMA 
best fit. In addition, as discussed in [031, the statistical analysis of KamLAND data shows no hints for 
anomalous effects beyond the standard scenario involving known reactor sources and neutrino flavor 
disappearance. 
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3.4 Solar and KamLAND constraints (#13 free) 

For #i3 > 0, electron neutrino mixing includes also u 3 (besides v\ and v 2 ), and Am 2 -driven oscillations 
can take place, with amplitude governed by #13. Therefore, the v e survival probability for # 13 > (P 3v ) 
generally differs from the one for 6*13 = (P 2v )- 

In KamLAND, Am 2 -driven oscillations are so fast to be smeared away by the finite energy resolution, 
leaving only the average (#13 mixing) effect, 

Pzv = 4 3 P 2v + 4 . (34) 

In first approximation, a similar formula holds for solar neutrinos, provided that the neutrino potential 
V is multiplied everywhere by cos 2 #13 (see |%T | I122j and refs. therein): 

Pzv * 4A + 4 , (35) 

Plv — ► P-iv = P2is\ V ^c( 3 V ■ fi6) 

This replacement generates a mild energy-dependence of the correction, which is absent in Eq. (J34j) . 

In second approximation, solar neutrinos develop a subleading dependence of P 3v on Am 2 and 
on its sign (i.e., on the hierarchy, see [HS]; such dependence disappears for Am 2 — > 00, where one 
recovers the above equations. Accurate analytic expressions for the subleading Am 2 effects on P 3v as a 
function of energy can be found in [SB] . 

Figure 13 shows the size of leading (6* 13 -driven) and subleading (±Am 2 -driven) effects, through the 
fractional difference between P 2v and P 3u , calculated for the representative value s 2 3 = 0.04 and for 
best-fit LMA parameters (and averaged over the 8 B solar neutrino production region, for definiteness) . 
The solid curve is calculated for Am 2 = 00, i.e., no subleading effect; the leading effect (about —7%) 
is almost entirely due to the factor cf 3 in front of P 2v , plus a mild energy dependence. The dashed and 
dot-dashed curves are instead calculated for Am 2 = +2.4 and —2.4 (xlO -3 eV 2 ), respectively; their 
difference from the solid curve quantifies the size of Am 2 subleading effects. Although the dependence 
of P 3v on Am 2 and on the hierarchy is theoretically interesting (see, e.g., |122j ). such subleading effect is 
an order of magnitude smaller than the "leading" ^13-effect in Fig. 13, and its inclusion would not change 
in any appreciable way the analysis of solar neutrino data (as we have explicitly checked). Therefore, 
in the following, we can safely assume the approximations in Eqs. (j33|) and JHSJ), i.e., neglect the effect 
of Am 2 and its sign in the solar(+KamLAND) data analysis, as it was the case for older data 123 . 

Figure 14 shows the results of our analysis of solar and KamLAND data (both separately and in 
combination) for unconstrained values of 9\ 3 , in terms of the 2a projections of the (5m 2 , sin 2 9% 2 , sin 2 #13) 
allowed region onto each of the three coordinate planes. There is no statistically significant preference 
for # 13 7^ 0, and upper bounds are placed by both solar and KamLAND data separately. 

Concerning KamLAND data only, there is a slight anticorrelation between sin 2 # 13 and sin 2 (9 12 (upper 
left panel in Fig. 14), since the total rate information constrains both parameters |124j . and a higher 
sin 2 #13 can be traded for a lower sin 2 9\ 2 . However, sin 2 Q\ 2 cannot decrease indefinitely — since it would 
suppress the amplitude of the observed shape distortions [20] — and thus an upper bound on sin 2 9\ 3 
emerges in KamLAND. 8 

We remind that the solar v sensitivity to sin 2 9\ 3 (Fig. 14) comes from the combination of all solar 
neutrino experiments, in contrast with the bounds on (5m 2 , sin 2 9i 2 ), which are dominated by the "high 
energy" 8 B neutrino experiments (SNO and SK). As discussed, e.g., in |122j . for increasing values of 
#13 a tension arises among different data sets and, in particular, between SNO and Gallium data. Such 
two experiments, probing respectively the high and low energy part of the solar neutrino spectrum, 

8 The KamLAND analysis in this work includes event-by-event energy information |21j but not the event time infor- 
mation |93j. We have explicitly checked that the time information, which does not significantly alter the bounds on 
(5m 2 , sin 2 #i 2 ) US], has also negligible effects on the the bounds on sin 2 6i 3 . 
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exhibit different correlation properties between the two mixing parameters 612 and #13. In particular, 
for increasing values of #13, the SNO and Gallium experiments tend to prefer higher and lower values 
of sin 2 # 12 , respectively |122| . worsening the good agreement currently reached at #13 ~ 0. Therefore, a 
"collective" effect of different experiments is responsible for the solar neutrino constraints on sin 2 #13. 
(See also [2U] for a discussion of bounds on 813 with earlier data.) 

Very interestingly, the combination of solar and KamLAND data in Fig. 14 is now powerful enough to 
place a combined upper bound on sin 2 #13 at the 5% level at 2a, not much weaker than the bound coming 
from the CHOOZ plus atmospheric data discussed below in Sec. 4.3 (see also |125j for an earlier dis- 
cussion of solar+KamLAND constraints on #13). Notice also that, in the combined (solar+KamLAND) 
regions of Fig. 14, there are negligible correlations among the three parameters; this fact implies that 
the bounds in Eqs. (|23|) and (|25j). derived for 13 = 0, hold without significant changes also for 13 
unconstrained. It also justifies (a posteriori) our choice to discuss in detail the case #13 = 0, which 
embeds most of the relevant information on the leading parameters (5m 2 , sin 2 # 12 ). 

4 SK atmospheric neutrinos, K2K, and CHOOZ 

In this Section we discuss the constraints on the mass-mixing parameters (Am 2 , 823, 613) coming from 
the SK atmospheric neutrino detector [26J, from the K2K long-baseline accelerator neutrino experiment 
[30, : 3TJ, and from the short-baseline CHOOZ reactor neutrino experiment [3*2*] . 

Our SK atmospheric neutrino analysis is performed by using the same event classification (binning) 
and systematic error treatment as in |117j . In particular, we consider (in order of increasing average 
energy) the zenith angle distributions of the so-called Sub-GeV (SG) electron and muon samples (SGe 
and SG//) in 10 zenith angle (8 Z ) bins; Multi-GeV (MG) electron and muon samples (MGe and MG//) in 
10 zenith angle (8 Z ) bins; Upward Stopping muons (US/i) in 5 bins; and Upward Through-going muons 
(UT/x) in 10 bins, for a total of 55 accurately computed observables. We include 11 sources of systematic 
errors |117j with the pull method (94 which allows a better understanding of systematic shifts. 9 With 
respect to our previous SK analysis |117j . we use updated results j2H] and — unless otherwise stated- 
atmospheric neutrino input fluxes from the three-dimensional (3D) simulation of 126] (see also |127H128] 
for other 3D results). 

Concerning the K2K experiment, we use the latest spectrum data from [ST], but regrouped in the 
same 6 bins as in |117j (by using information from |129j ); this choice is motivated by the fact that 
information about K2K correlated systematics has been made publicly available only for 6 bins (see 
|117j and references therein). Finally, the CHOOZ spectral data [32] are analyzed as in [S3]. Further 
technical details are given in the Appendix. 

In the following, we discuss first the impact of the SK+K2K data on the neutrino parameters 
(Am 2 , sin 2 #23) for 6*13 = 0. This allows to appreciate the sub leading effect induced by nonzero values 
of (5m 2 , sin 2 #12), especially on atmospheric neutrinos. Then we consider the more general case 6x3 7^ 0, 
and discuss in some detail the related subleading effects in SK, as well as the constraints from the 
SK+K2K+CHOOZ analysis. 

A final remark is in order. The MACRO [27] and Soudan-2 [28J atmospheric neutrino experiments 
provide (Am 2 , sin 2 #23) constraints which are consistent with those from SK [2B], but are also affected by 
larger uncertainties (due to the lower statistics and narrower L/E range); they are not included in this 
work. Similarly, the negative results of the Palo Verde reactor experiment [33] and of the K2K searches 
in — > u e appearance mode (.130] (consistent with, but less constraining than CHOOZ 32 ) are not 
included here. Future improved global analyses might take into account these additional data, finer SK 

9 The SK Collaboration has used a finer classification of events and systematics in [22], as well as an alternative L/E 
binning in [23]. Such refined analyses cannot be performed outside the Collaboration. 
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and K2K spectral binning, and the covariance of the SK and K2K common systematics (interaction 
cross section, detector fiducial volume, and event reconstruction errors). 

4.1 SK and K2K constraints for #13 = and statistical checks 

While for # 13 = the solar+KamLAND parameter space reduces to (5m 2 , sin 2 # 12 ) exactly, the atmo- 
spheric+K2K parameter space reduces to (Am 2 , sin 2 6*23) only to a first approximation. Indeed, while 
the assumption # 13 = forbids solar and reactor v e — > v^ T transitions involving u 3 and its associated 
parameters (Am 2 , sin 2 #23) , it does not forbid, e.g., atmospheric z/^ — > v T transitions involving the pair 
(z/1,^2), which depend on the 5m 2 parameter. This is most easily seen in the vacuum case, where 
Eq. (j2*U|) implies that, for i,j = 1,2 and even for 813 = 0, the P MT transition probability contains the 
following nonzero (v\, i^-induced term, 

- W, xUta U T ,U T2 sin 2 -44444 sin 2 l) . (37) 

The small effect of nonzero (LMA) values of (5m 2 , sin 2 #12) in the atmospheric neutrino data analysis, 
phenomenologically noted in [T23J for any 6*13, has often been legitimately neglected (except occasionally, 
see the bibliography in |131j ). being basically hidden by large statistical and systematic uncertainties. 
The full implementation of such effect is nontrivial (it requires a numerical 3z/ evolution in the Earth 
matter layers), and its main theoretical aspects have been elucidated only recently |131| 1132] . in con- 
nection with the progressive confirmation and determination of the LMA parameters by solar and 
KamLAND data, and with the increasing accuracy of atmospheric neutrino data. Although still small, 
the effect is definitely not smaller than others which are usually taken care of, and deserves to be 
included in state-of-the-art analyses |132l 11331 H34j . 

For instance, Fig. 15 shows the results of our analysis of the latest SK data in the plane (Am 2 , sin 2 #23) 
at #13 = 0, for three increasingly accurate inputs: Atmospheric neutrino fluxes from one-dimensional 
(ID) simulations and 5m 2 = (top panel); atmospheric neutrino fluxes from full three-dimensional 
(3D) simulations [ 126| and 5m 2 = (middle panel); and finally, 3D fluxes and LMA best-fit values 
for (5m 2 , sin 2 12 ). In all panels, the three curves refers to 1, 2, and 3a contours, and the best-fit 
point is marked by horizontal and vertical lines to guide the eye. One can appreciate that the (now 
customarily included) 3D flux input shifts Am 2 downward by ~ 0.5a with respect to the ID flux input; 
on the other hand, the inclusion of subleading LMA effects shifts sin 2 623 by ~ 0.5<r |133j with respect 
to the hypothetical case 5m 2 = 0. As expected, both effects are small, but there is no reason to 
keep the first and to neglect the second. Moreover, the LMA effect intriguingly breaks the 823 octant 
sysmmetry, which is in principle an important indication for model building (see, e.g., |135j ). Hereafter, 
the analysis of the SK atmospheric data will explicitly include nonzero values of (5m 2 , sin 2 #12), fixed at 
their best-fit values in Eqs. (|25 | 12 6 |) but with no uncertainty (whose effect is really negligible). For the 
sake of completeness, LMA-induced and matter effects will also be included in the calculation of the 
K2K oscillation probabilities, where, however, such effects are even smaller than in SK, as discussed in 
Sec. 4.2. 

Figure 16 shows, for 9 13 = 0, the results of our analysis of SK and K2K data, both separately and in 
combination. Notice that the top panel in Fig. 16 is the same as the bottom panel in Fig. 15. The K2K 
constraints are octant-symmetric and relatively weak in sin 2 623, while they contribute appreciably to 
reduce the overall Am 2 uncertainty. Therefore, not only K2K confirms the neutrino oscillation solution 
to the atmospheric neutrino anomaly with accelerator neutrinos j3T|, but it also helps in reducing the 
oscillation parameter space. Moreover, there is still room for improvements in K2K. Figure 17 shows 
our both unoscillated and oscillated K2K spectrum of events (at the SK+K2K best-fit in Fig. 16) in 
terms of the reconstructed neutrino energy, as used in this work. The oscillated spectrum is shown both 
with and without the systematic shifts in our pull approach; such shifts are modest as compared with 
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the large statistical errors. Therefore, one can reasonably expect that, with higher statistics, the final 
K2K data sample can further contribute to reduce the Am 2 uncertainty. 

Systematic effects are instead quite important in the SK atmospheric neutrino analysis. Fig. 18 shows 
the ratio of experimental data 26J and of best-fit theoretical predictions (with and without systematic 
pulls) with respect to no oscillations, as a function of the zenith angle of the scattered lepton (e or /i), 
for the five samples used in the analysis. In particular, in terms of Eq. (}3*3*j) . the dashed histograms 
represent the unshifted theoretical predictions (central values Rn), while the dashed histograms represent 
the systematically shifted predictions (R n ) for the given mass-mixing parameters (which correspond to 
the SK+K2K best-fit point in Fig. 16). Vertical error bars represent the la statistical uncertainties of 
the data. The electron data (SGe and MGe) show some excess with respect to the unshifted predictions 
(dashed lines), which tends to be reduced when systematic shifts are allowed (solid lines). Notice that 
the dashed lines slightly differ from unity for upward (cos 6^ ~ —1) events in the SGe and MGe sample, 
as a result of subleading LMA effects (5m 2 ^ 0). A systematic, upward shift of the predictions is also 
preferred in the high-energy muon samples, US/i and UT/i, and especially in the latter, where it amounts 
to ~ 20%. The pull analysis of the observables in Fig. 19 tells us that such shifts are not necessarily 
alarming from a statistical viewpoint, since they are all smaller than two standard deviations. However, 
their distribution is definitely not random: within each of the six data samples, most of the pulls in 
Fig. 19 are one-sided, indicating that there seems to be some normalization offset. This is confirmed by 
the pull analysis of the systematics in Fig. 20, where the two largest pulls (~ 1.5<r) refer to normalization 
parameters (p and pt) which govern the relative normalization of muon samples with increasing energy 
(fully-contained, partially-contained, and upward-stopping muons) [1 17j . Also the sub-GeV muon-to- 
electron flavor ratio error (/3 S ) is stretched beyond la in Fig. 20. Although there is no alarming "3<r" 
offset anywhere, it is clear that a better understanding and reduction of the systematic error sources 
(i.e., atmospheric neutrino fluxes, interaction cross sections, detector uncertainties) is needed (134] if 
one wants to observe in the future small subleading effects, as those induced by Sm 2 ^ and #13 7^ 
and discussed in more detail in the next section. 

4.2 Discussion of subleading effects 

Our calculations of atmospheric neutrino oscillations are based on a full three-flavor numerical evolution 
of the Hamiltonian along the neutrino path in the atmosphere and (below horizon) in the known Earth 
layers |123t I136| 11371 H38j . Semianalytical approximations to the full numerical evolution (although 
not used in the final results) can, however, be useful to understand the behavior of the oscillation 
probability and of some atmospheric neutrino observables. A particularly important observable is the 
excess of expected electron events (N e ) as compared to no oscillations (iV°): 



where P a p = P(v a — > z/^), and r is the ratio of atmospheric and v e fluxes (r ~ 2 and ~ 3.5 at sub-GeV 
and multi-GeV energies, respectively). In fact, this quantity is zero when both # 13 = and 5m 2 = 0, 
and is thus well suited to study the associated subleading effects (which may carry a dependence on 
the matter density) in cases when Sm 2 and #13 are different from zero |131j . 

We remind that matter effects are governed by A(x) = 2^/2GfN c (x)E, with N e ~ 2 mol/cm 3 in the 
Earth mantle and ~ 5 mol/cm 3 in the core. It can be easily derived that 



implying that Earth matter can substantially affect Am 2 -driven oscillations [i.e., A/ Am 2 ~ 0(1)] for 



l = (P ee -l)+rP t 



(38) 




(39) 
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E ~ 0(10) GeV, i.e., in multi-GeV and upward-stopping events. Similarly, 
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3.8 



8 x 10" 



5m 2 




E 



1 GeV J \2 mol/cni 



(40) 



ie and P e/I can be 



implying that A /8m 2 ~ 0(1) for sub-GeV SK events (and, in principle, for accelerator K2K neutrinos as 
well). In the constant-density approximation A(x) = const (i.e., by neglecting mantle-core interference 
effects |138j to keep the following discussion simple), the oscillation probabilities P ( 
evaluated through Eq. (|21|). in terms of the effective mass- mixing parameters in matter ^ v%3 , 

Suitable approximations for such parameters have been reported in many papers. If we use, e.g. 
those reported in the classic review [HS], after some algebra we get from Eqs. (|HK|) and ()21)1 that the 
electron excess at sub- or multi-GeV energies can be written as a sum of three terms, 
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(46) 



The above expressions for Aj, which hold for neutrinos with normal hierarchy and 5 = 0, coincide 
with those reported in |131j (up to higher-order terms or CP- violating terms, not included here). The 
corresponding expressions for antineutrinos, for inverted hierarchy, and for 5 = it, can be obtained, 
respectively, through the replacements: 



+ A -> -A 
-Am 2 — > —Am 2 

+Sl3 ~~ > ~ s 13 



(swaps (anti)neutrinos) , 
(swaps hierarchy) , 
(swaps CP parity) , 



(47) 
(48) 
(49) 



where by "CP parity" we mean cos 5 = ±1. Under such transformations, the terms Aj behave as follows: 

(1) all Aj's are affected by A —>■ —A through 12 or 13 ; (2) only A x is sensitive to Am 2 — ► —Am 2 ; (3) 
only A 3 is sensitive to +Si 3 — > — s 13 . 

Concerning the dependence on the oscillation parameters, one has that: (1) all Aj's depend on 2 3; 

(2) Ai arises for 13 > 0, and is independent of 5m 2 ; (3) A 2 arises for 5m 2 > 0, and is independent of 
013 ; only A 3 ("interference term" |131p depends on both i3 and 5m 2 . 

Concerning the dependence on energy, in the sub-GeV range one has that: (1) 0i3 ~ 0i 3 , so that 
for large L the first term is simply Ai ~ 2sf 3 Ci 3 (rs2 3 — 1); (2) since r ~ 2, the term Ai flips sign as 
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s| 3 crosses the maximal mixing value 1/2 |139j . and similarly for A2 (with opposite sign) |131j : (3) for 
neutrinos, which give the largest contribution to atmospheric events, it turns out that tan2# 12 < 0, and 
thus typically A 3 < for 5 = (A 3 > for 5 = ir). In the multi-GeV range one has that 9 U — 7r/2, 
so that only Ai dominates, with typically positive values (being r ~ 3.5 and s 23 not too different from 
1/2). 

Figure 21 shows exact numerical examples (extracted from our SK data analysis) where, from top 
to bottom, the dominant term is Ai, A2, and A3. Here, as in Fig. 18, the dashed histograms represent 
the unshifted theoretical predictions, while the solid histograms represent the systematically shifted 
predictions, i.e., R n and R n respectively [in terms of Eq. (JHBJ)]- Let us focus on subleading effects in the 
dashed histograms of Fig. 21, which refer to the sub-GeV (left) and multi-GeV (right) electron samples. 
In the figure, we have taken Am 2 = +2.4 x xlO -3 eV 2 (normal hierarchy); other relevant parameters 
are indicated at the right of each panel. In the upper panel, we have set 5m 2 = 0, so as to switch off 
A 2 and A 3 . We have also taken s| 3 = 0.4 < 0.5, so that A x < in the sub-GeV sample; it is instead 
Ai > in the multi-GeV sample. In the middle panel, we have set (5m 2 , sin 2 # 12 ) at their best-fit LMA 
values, but have taken sin 2 #13 = 0, so that only A 2 survives. In particular, while there is no observable 
effect of A 2 in the multi-GeV sample (where the energy is relatively high and sin2#i 2 ~ 0), the effect is 
positive for sub-GeV neutrinos, where s| 3 = 0.4 < 1/2. Notice that the upper and middle panel results 
are insensitive to 5 = or tt, since A 3 ~ in both cases. Finally, in the bottom plot we have taken 
s 23 = 1/2, so as to suppress Ai and A 2 is the sub-GeV sample, where A 3 > for our choice 5 = tt. In 
the multi-GeV sample, however, Ai is still operative. 

The subleading dependence of atmospheric electron neutrino events on the hierarchy, 5m 2 , 6*i 3 , and 
CP-parity is intriguing and is thus attracting increasing interest |140j . However, Fig. 21 clearly shows 
that such dependence is currently well hidden, not only by statistical uncertainties (vertical error bars) 
but, more dangerously, by allowed systematic shifts of the theoretical predictions (solid histograms). 
For instance, in the upper panel, systematics can "undo" the negative effect of Ai in the SGe sample 
and make it appear positive. In all cases, they tend to magnify the zenith spectrum distortion; this is 
particularly evident in the right middle panel, where the unshifted theoretical prediction is flat. 

We think it useful to quantify at which level one has to reduce systematic uncertainties, in order to 
appreciate subleading effects in future, larger SK-like atmospheric neutrino experiments such as those 
proposed in |141| 1142] 1143] (see also |144| ). Since normalization systematics are large (as discussed in 
the previous section) and a significant reduction may be difficult, we prefer to focus on a normalization- 
independent quantity, namely, the fractional deviation of the up-down asymmetry of electron events 
from their no-oscillation value, 

where "up" (U) and "down" (D) refer to the zenith angle ranges cos #2 G [—1,-0.4] and [0.4,1], 
respectively. We perform a full numerical calculation of this quantity for both SGe and MGe events, 
assuming the SK experimental setting for definiteness. Notice that the up-down asymmetry involves 
the first and last three bins of the SGe and MGe samples in Fig. 18. 

Fig. 22 shows isolines of 100 x A e for the SGe sample, plotted in the (sin 2 # 23 , sin 2 # 13 ) plane at fixed 
Am 2 = 2.4 x 10~ 3 eV 2 , for both normal hierarchy (+Am 2 , left panels) and inverse hierarchy (—Am 2 , 
right panels). In both hierarchies, we consider first the "academic" case 5m 2 = (top panels), then we 
switch on the LMA parameters (5m 2 , sin 2 # 12 ) at their best-fit values, for the the two CP-conserving 
cases 5 = (middle panels) and 5 = tt (bottom panels). The isolines in the upper panels reflect the 
behavior of the Ai term, which is positive (negative) for s 23 > 1/r (< 1/r), and vanishes for #13 — > 0, 
with a weak dependence on the hierarchy through #13. In the middle panel, subleading LMA effects are 
operative through A 2 and A3. The variation of A e in sign and magnitude is now more modest as s| 3 
increases, since the variation of the Ai term is now partially compensated by the opposite variation of 
the A 2 term. In particular, the term A 2 is responsible for nonzero values of A e at 9 is = 0, which break 
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the 823 octant symmetry |133j . as also phenomenologically observed in Fig. 15. The difference between 
the middle and bottom panels is due to the interference term A 3 , which is typically negative (positive) 
for 5 = (5 = 7r), and thus either adds or subtracts to A x and A 2 in the two cases. However, for # 13 = 
the values of A e basically coincide in all middle and bottom panels, since the only surviving term (A 2 ) 
carries no dependence on the hierarchy or the CP parity. From the results in Fig. 22 we learn that: 
(1) subleading <5m 2 -induced effects are of the same size of ^-induced effects in the SG sample, so none 
can be neglected in a precise 3z/ oscillation analysis; (2) for nonzero values of both 5m 2 and #13, the 
sub-GeV electron asymmetry is typically more pronounced (and positive) for 5 = it, as compared with 
the case 5 = 0; one can thus expect the latter case to be slightly disfavored in a global fit (since the 
SGe data show a slight asymmetry, see Fig. 18); (3) in any case, the electron asymmetry is typically at 
the percent or sub-percent level for sin 2 # 13 < few%; therefore, statistical and systematic uncertainties 
need to be reduced at this extraordinary small level in order to really "observe" the effects in future 
atmospheric neutrino experiments [134 . 

Fig. 23 shows our numerical calculation of the up-down electron asymmetry for the SK multi-GeV 
sample. The six panels refer to the same cases as in Fig. 22. In the MGe sample, the terms A 2 and 
A3 are small, and there is little dependence on 5m 2 and on the CP parity (top, middle and bottom 
panels being quite similar). The dominant Ai term makes the asymmetry generally positive, and with 
significant dependence on the hierarchy (left vs right panels) through # 13 . The MGe asymmetry can 
be of 0(10%) and thus relatively large; with some luck, such asymmetry might be seen in future large 
Cherenkov detectors if #13 is not too small (see, e.g., |145j and refs. therein). In any case, one can expect 
some dependence of the current SK fit on the hierarchy through multi-GeV events (see also |137j for 
older data); it is difficult, however, to "predict" which of the two hierarchies (normal or inverted) is 
currently preferred, since large statistical fluctuations make the zenith-angle pattern of SK MGe data 
somewhat erratical (see Fig. 18). 

We conclude this Section with a brief note on subleading effects in K2K (which have been numerically 
included throughout this work). For #13 = 0, it is easy to derive that, in vacuum, 

9 ~ f Am 2 - ^fcos6 12 \ , . 

p w — 1 - sin 2#23 sin f ^ \ , (51) 

at first order in 5m 2 / Am 2 . It is also not difficult to check that this formula is not significantly affected 
by matter effects, as far as A > 5m 2 , which is true for most of the K2K event spectrum. The #23 
octant-symmetry of the above equation is responsible for the appearance of two degenerate best-fits in 
the K2K analysis of Fig. 16. The above equation is not invariant under a change of hierarchy [146 
(+Am 2 — > —Am 2 ), which leads to a (really tiny) relative change of the oscillation phase equal to 
±(5m 2 /2) cos 6 12/ Am 2 ~ 0.6%; this change does not produce graphically observable effects in Fig. 16. 
Although very small, these and other subleading K2K effects (e.g., those arising for both #13 and 
5m 2 nonzero) have been kept in the analysis, in order to be consistent with the atmospheric neutrino 
data analysis (where such effects have also been included, as already discussed), and in order to show 
explicitly their impact on the global SK+K2K+CHOOZ analysis presented in the next section. 



4.3 SK, K2K and CHOOZ constraints (0 13 free) 

In this section we present the results of our analysis of SK+K2K+CHOOZ data for unconstrained values 
of (Am 2 , sin 2 #23, sin 2 6 13 ), and for fixed values (5m 2 , sin 2 6> 12 ) = (8 x 10~ 5 eV 2 , 0.314) in all the three 
data samples. There are four discrete subcases in our analysis, corresponding to a change in hierarchy 
or CP parity: 

[sign(Am 2 ) = ±1] <g> [cos 5 = ±1] . (52) 
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In particular we remind that, in this work, we do not consider generic values of 5, but only the 
two inequivalent CP-conserving cases (5 = and 8 = ir). They are related by the transformation 
+ s i3 - ► _ Si3 which, of course, does not mean that s 13 can be negative, but just that cos5si 3 can 
change sign. Since the two cases smoothly merge for S13 — > 0, we think it useful to show the results 
of our analysis also in terms of the variable cos 5 sin #13, i.e., of ± sin #13, for both normal and inverted 
hierarchy. 

Figure 24 shows the x 2 function from the SK+K2K+CHOOZ fit, in terms of cos<5 sin #13, for 
marginalized (Am 2 , sin 2 # 23 ) parameters. 10 The solid and dashed curves correspond to normal and 
inverted hierarchy, respectively, while their left and right parts correspond to 5 = n and 5 = 0, respec- 
tively. Notice that the solid and dashed curves do not exactly coincide at 6*13 = 0, since for 5m 2 > 
there is a very weak dependence on the hierarchy even at #13 = in reactor [551 1147] . accelerator .146J, 
and atmospheric |148j neutrino oscillations. The difference is, however, really tiny within the current 
global analysis (A% 2 ~ 0.2 at s 13 = 0). The absolute x 2 minimum is reached in the left half of the figure 
(5 = it) for sin #13 ~ 0.1; the minimum in the right half (5 = 0), which is reached for #13 = 0, is only 
slight higher (A% 2 < 1). The slight difference between these two CP-conserving cases is mainly due to 
sub-GeV SK events, as discussed in the comments to Fig. 22. Finally, normal and inverted hierarchies 
give basically the same results for small values of S13 (say, < 0.1), while the latter hierarchy is slightly 
preferred for higher values of S13. The fit becomes rapidly worse for S13 ~ 0.2 or higher. Figure 24 
nicely summarizes our current (unfortunately weak) sensitivity to the neutrino mass hierarchy and to 
the extremal (CP-conserving) cases 5 = and 5 = 71. 

Figure 25 shows the parameter space orthogonal to the one in Fig. 24, i.e., the bounds on (Am 2 , sin 2 #23) 
for marginalized #13, in each of the four cases in Eq. (|52j) . The differences between such cases are very 
small. Figure 24 and 25 confirm that our current sensitivity to the subleading effects — which distinguish 
the four subcases in Eq. (J52*Jl — is not statistically appreciable yet. Therefore, it makes sense to make 
a further marginalization over the four subcases, by minimizing the SK+K2K+CHOOZ function with 
respect to hierarchy and CP parity. The results are shown in Fig. 26, in terms of the projections of 
the (Am 2 , sin 2 6*23, sin 2 6*13) region allowed at 1, 2, and 3a onto each of the coordinate planes. The best 
fit is reached for nonzero 9 13 (as expected from Fig. 24), but #13 = is allowed within less than la. 
The preferred value of sin 2 623 remains slightly below maximal mixing. The best-fit value of Am 2 is 
2.4 x 10 -3 eV 2 . Notice that the correlations among the three parameters in Fig. 26 are very weak. 

5 Global analysis of oscillation data 

The results of the global analysis of solar and KamLAND data (Sec. 3.4) and of SK+K2K+CHOOZ data 
(Sec. 4.3) can now be merged to provide our best estimates of the five parameters (5m 2 , Am 2 , 12 , #13, 623), 
marginalized over the four cases in Eq. (|52p. The bounds will be directly shown in terms of the "number 
of sigmas", corresponding to the function (A% 2 ) 1//2 for each parameter. 

Figure 27 shows our global bounds on sin 2 613, as coming from all data (solid line) and from the 
following partial data sets: KamLAND (dotted), solar (dot-dashed), solar+KamLAND (short-dashed) 
and SK+K2K+CHOOZ (long-dashed). Only the latter set, as observed before, gives a weak indication 
for nonzero 613. Interestingly, solar +KamL AND data are now sufficiently accurate to provide bounds 
which are not much weaker than the dominant SK+K2K+CHOOZ ones, also because the latter slightly 
prefer 613 > as best fit, while the former do not. 

Figure 28 shows our global bounds on the four mass-mixing parameters which present both upper 
and lower limits with high statistical significance. Notice that the accuracy of the parameter estimate 
is already good enough to lead to almost "linear" errors, especially for 5m 2 and sin 2 # 12 . For Am 2 and 

10 This representation is inspired by Ref. |132j where, however, cos<5sm 2 # 13 was used. The use of cos<5sin#x3 makes 
the x 2 curves smooth (no cusp) across #13 = 0. 
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sin 2 #23, such "linearity" is somewhat worse in the region close to the best fit (say, within ±ler), and 
thus 2a (or 3a) errors should be taken as reference. 

We summarize our results through the following ±2cx ranges (95% CL.) for each parameter: 



Notice that the lower uncertainty on sin 2 #13 is purely formal, corresponding to the positivity con- 
straint sin 2 #13 > 0. Correlations among parameters are not quoted, being currently small (as already 
observed) . 

The above bounds have been obtained from a global analysis of oscillation data (for U = U*). 
They have, however, an impact also on non-oscillation observables. In particular, the smallness of the 
squared mass splittings induces significant correlations on the three parameters (mp,mpp, E) which are 
sensitive to absolute neutrino masses (see [S3] and references therein). Figure 29 shows the updated 
2a allowed bands in each of the three corresponding coordinate planes, for both normal and inverted 
hierarchy. There is an evident positive correlation, especially between mp and E; the correlation is less 
pronounced when it involves mpp, due to our ignorance of the Majorana phases 02 and 03 (that we 
take as free parameters). The two hierarchies split up only at very low values of the observables, where 
mass splittings start to be of the order of the absolute masses (non-degenerate cases). Non-oscillation 
data on mp, mpp and E can reduce the allowed parameter space in Fig. 29, hopefully leading to a single 
solution and thus to the determination of the absolute neutrino masses. Such data are discussed in the 
next Section. 

6 Global analysis of oscillation and non-oscillation data 

In this Section we discuss first non-oscillation data on the three observables {mp, mpp, E), and then show 
how these data further constrain and reduce the allowed regions in Fig. 29. As we shall see, when all the 
data are taken at face value, no combination is possible: a strong tension arises, indicating that either 
some experimental information or their theoretical interpretation is wrong or biased. In particular, it 
appears difficult to reconcile [53] the claimed Oz/2/3 signal [17] and the most recent upper bounds on E 
from precision cosmology |4*T] I42[ 14*3*] . However, relaxing one of either pieces of data reduces the tension 
and allows a global combination, which can be valuable for prospective studies [55]. 

Needless to say, the relations between the variables (mp,mpp,Tj) have been subject to intensive 
studies, which form a large specialized literature on absolute neutrino mass observables. We refer the 
reader to the review papers jSHl HUB ESI E7J EHl I149j for extensive bibliographies, and to the articles 
jSHl E7J H5Uj for recent up-to-date discussions. 

6.1 Bounds on mp 

Experimental constraints on the effective electron neutrino mass mp have been recently presented [3B] 
for the Mainz and Troitsk tritium /3-decay experiments. The experimental values are consistent with 
zero within errors. Their combined upper bound at 2a has been estimated in [55] as: 



which is less conservative than the 3 eV upper limit recommended in pQ. It should be mentioned 
that the Troitsk results are to be taken with some caution, being affected by an unexplained anomaly 




0.9±H x 10" 2 , 
7.92 (1 ± 0.09) x 10~ 5 eV 2 
0.314 (li^) , 
2.4 (l+_° il) x 10- 3 eV 2 , 
0.44 . 



(53) 
(54) 
(55) 
(56) 
(57) 



mp < 1.8 eV (Mainz + Troitsk) , 



(58) 
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(namely, a fluctuating excess of counts near the endpoint) [SH]. However, as we will see, upper limits 
on mp in the 2-3 eV range are, in any case, too weak to contribute significantly to the current global 
fit in the (mp, mpp, S) parameter space, so that "conservativeness" is not (yet) an issue in this context. 



6.2 Bounds on mpp 

Neutrinoless double beta decay processes of the kind (Z, A) — > (Z + 2, A) + 2e~ have been searched in 
many experiments with different isotopes, yielding negative results (see JHHIC49J for reviews). Recently, 
members of the Heidelberg-Moscow experiment have claimed the detection of a Oz/2/3 signal from the 
76 Ge isotope Wf\. If this signal is entirely due to light Majorana neutrino masses, the Oz/2/3 half-life 
T is related to the mpp parameter by the relation 

2 

m%p = , (59) 

'—'mm-' 

where m e is the electron mass and C mm is the nuclear matrix element for the considered isotope [3*H] . 

Unfortunately, theoretical uncertainties on C mm are rather large (see e.g. [3"9*]). and their — somewhat 
arbitrary — estimate is matter of debate (see |151|. 11521 1153j and refs. therein). In [S3] we adopted a 
naive but very conservative estimate, by defining the range spanned by "extremal" published values of 
Cmm as an "effective 3a range," thus obtaining log 10 (C mm /y~ 1 ) = —13.36 ± 0.97 (at ±3cr). Here we 
prefer to adopt the results of a recent detailed discussion of the nuclear model uncertainties for C mm , 
performed within the (Renormalized) Quasiparticle Random Phase Approximation, and calibrated to 
known 2v(5(5 decay rates |154j . For our purposes, we cast the results of such promising approach |154j in 
the form log 10 (C , mm /y _1 ) = —13.36 ± 0.15 (at ±3cr), where systematic coupling constant uncertainties 
(gA = 1-1.25, see |154j ) have been included. This "new" range for C mm has (accidentally) the same 
central value as before, but with significantly reduced errors. Under the assumption of a positive 0z/2/5 
signal jUj, we then derive that 

log 10 (m^/eV) = -0.23 ± 0.14 (2a) , (60) 

i.e., 0.43 < mpp < 0.81 (at 2a, in eV). See also jSS] for our previous (more conservative) estimated 
range. 

The claim in H7j has been subject to strong criticism, especially after the first publication [33] 
(see |391 1149j and refs. therein). Therefore, we will also consider the possibility that T = oo is allowed 
(i.e., that there is no 0^2/3 signal), in which case the experimental lower bound on mpp disappears, and 
only the upper bound remains. In conclusion, we adopt the following two possible Oz/2/? inputs for our 
global analysis: 

log 10 (m w /eV) = -0.23 ± 0.14 (0u2(3 signal assumed) , (61) 
log 10 (m w /eV) = -0.23±^ 14 (0u2/3 signal not assumed) , (62) 

where errors are at 2a level. Concerning the unknown Majorana phases 02 and 03 in Eq. (J15)) . we 
simply assume that they are independent and uniformly distributed in the range [0, it], which covers all 
physically different cases in mpp. 



6.3 Bounds on E 

The neutrino contribution to the overall energy density of the universe can play a relevant role in large 
scale structure formation, leaving key signatures in several cosmological data sets. More specifically, 
neutrinos suppress the growth of fluctuations on scales below the horizon when they become non rel- 
ativists. Massive neutrinos of a fraction of eV would therefore produce a significant suppression in 
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the clustering on small cosmological scales. Data on large scale structures, combined with Cosmic 
Microwave Background (CMB) and other precision astrophysical data, can thus constrain the sum of 
neutrino masses E (see jlU l%9~l ITB31 11561 H57| for recent reviews). 11 

In this work we use the bounds on E previously obtained in collaboration with other authors in 
[55J, to which we refer the reader for technical details. We briefly remind that the experimental input 
used in |55J included CMB data from the Wilkinson Microwave Anisotropy Probe (WMAP) |41j . large 
scale structure data (221 from the 2 degrees Fiels (2dF) Galaxy Redshift Survey |158j and, optionally, 
constraints on mall scales from the recent Lyman a (Lya) forest data of the Sloan Digital Sky Survey 
(SDSS) |159| . The latter data have a strong impact on the current upper bounds on E (33], but are 
also affected by large systematics, which deserve further study |4lH 1159j . As in j^S], we conservatively 
quote (and use) upper bounds on E both with and without such Lya forest data; in particular, the 2a 
upper bounds from [S3] read: 



6.4 Impact of non-oscillation observables 

The experimental limits on the non-oscillation observables (rap, rapp, E), previously reported in terms 
of 2a ranges, are appropriately combined with oscillation data through A% 2 functions Although 
such combination can provide allowed regions at any confidence level, in the following we shall continue 
to show only 2a bounds for simplicity. 

Figure 30 shows the impact of all the available non-oscillation data, taken at face value, in the 
parameter space (rapp,Ti). Bounds on the third observable rap are projected away, being too weak 
to alter the discussion of the results in this figure. The horizontal band is allowed by the positive 
0z/2/5 experimental claim j32| equipped with the nuclear uncertainties of |154j through Eq. (|fil|). The 
slanted bands (for normal and inverted hierarchy) are allowed by all other neutrino data, i.e., by the 
combination of neutrino oscillation constraints (as in Fig. 29) and of astrophysical CMB+2dF+Lya 
constraints through Eq. ([63)1 . The tight cosmological upper bound on E prevents the overlap between 
the slanted and horizontal bands at 2a, indicating that no global combination of oscillation and non- 
oscillation data is possible in the sub-eV range. The "discrepancy" is even stronger than it was found 
in (53], due to the adoption of smaller 0^2/3 nuclear uncertainties [154J. It is premature, however, to 
derive any definite conclusion as to which piece of the data or of the 3u scenario is "wrong" in this 
conflicting picture. Further experimental and theoretical research is needed to clarify absolute neutrino 
observables in the sub-eV range. 

It is tempting, however, to see if the removal of some pieces of data can relax the tension in Fig. 30. 
The effect of removing only the lower bound on rapp through Eq. ([6"2~)) is shown in Fig. 31. Of the three 
remaining upper bounds on rap, on mpp, and on E, the latter is definitely dominant, and implies that 
future beta and double-beta decay searches should push their sensitivity below 0.2 eV, irrespective of 
the hierarchy. Conversely, the effect of removing only the Lya forest data through Eq. ([M)) is shown 
in Fig. 32. In this case, the combination of the claimed 0z/2/5 signal with oscillation data dominates 
the global fit, and "predicts" the observation of E ~ 1.5 eV and mp ~ 0.5 eV within formally small 
uncertainties (about ±20% at 2a). These predictions would really be "around the corner" from the 
observational viewpoint, both for E jBH] and for mp Future searches are expected to clarify the— 
currently controversial — situation about absolute mass observables in the sub-eV range, as depicted in 



11 Future cosmological data might become slightly sensitive to finer details (e.g., the neutrino mass hierarchy) through 
subleading effects |55] not included in this work. 



E < 0.5 eV (with Lya data) , 
E < 1.4 eV (without Lya data) . 



(63) 
(64) 



Figs. 30-32. 
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7 Conclusions 

We have performed a comprehensive phenomenological analysis of a vast amount of data from neutrino 
flavor oscillation and non-oscillation searches, including solar, atmospheric, reactor, accelerator, beta- 
decay, double-beta decay, and precision astrophysical observations. In the analysis, performed within 
the standard scenario with three massive and mixed neutrinos (for both mass hierarchies and for the 
two inequivalent CP-conserving cases), we have paid particular attention to implement subleading 
oscillation effects in numerical calculations, and to carefully include all known sources of uncertainties 
in the statistical comparison of theoretical predictions and experimental data. We have discussed 
the impact of solar and reactor data on the parameters (5m 2 , sin 2 #12, sin 2 #13), as well as the impact 
of atmospheric and reactor data on (Am 2 , sin 2 #23, sin 2 #13) . The bounds from the global analysis of 
oscillation data have been summarized, and several subleading effects have been discussed. Finally, 
we have analyzed the interplay between the oscillation parameters (5m 2 , Am 2 , sin 2 #12, sin 2 #23, sin 2 #13) 
and the non-oscillation observables sensitive to absolute neutrino masses (m^, m^, E), both with and 
without controversial data, which may or may not allow a reasonable global combination of all data. 
The detailed results discussed in this review article represent a state-of-the-art, accurate and up-to- 
date (as of August 2005) overview of the neutrino mass and mixing parameters within the standard 
three-generation framework. 
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Appendix 

In this Appendix we report — for more expert readers — additional technical information about the \ 2 
analysis of each data sample, which has been used to derive partial and global parameter bounds in 
this work. In particular, global constraints have been obtained by adding up all \ 2 contributions and 
by scanning the (CP-conserving) 3z/ mass-mixing parameter space 



with allowed regions being derived through A% 2 cuts with respect to the best-fit point. As discussed in 
the preceding Sections of this review, some data samples are actually not sensitive to all of the above 
parameters; the relevant variables in p will be explicitly emphasized in the following subsections. 



The CHOOZ experiment [32] has measured the positron energy spectra induced by 77 e 's produced by 
two nuclear reactors located at L\ = 1114.6 km and L 2 = 997.9 km from the detector. Each of the 
two spectra is divided in 7 energy bins, for a total of 14 event rate bins. In our analysis, such data are 
included through the following x 2 function: 



p = {±Am 2 , 5m 2 , s 2 23 , s 2 13 , s\ 2 , cos 5 = ±1} 



(65) 



7.1 CHOOZ 




,expt 




,cxpt 




(66) 
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where R^ xp and R\ are the observed and predicted rates in each bin, respectively, and a is an overall 
normalization factor with uncertainty a a = 2.7 x 10 -2 . The squared error matrix is defined as [3*2] : 

<4 = hi s 1 + u 1) + (kj-7 + kj+7)c 2 i , (67) 

where and Ui represent statistical errors and uncorrelated systematic errors, respectively, while the 
q's represent fully correlated systematic errors between equal-energy bins in the two reactor spectra. 
The theoretical rate in each bin is estimated as 

Rf e ° = ■ (Pee(p))i , (68) 

where R® is the no-oscillation rate, and (P ee (p))j is the V e survival probability, averaged over the 
appropriate energy range for the i-th bin, taking into account the detector energy resolution and the 
reactor core size. Numerical values for the rates R® and the errors Sj, Ui, and q, can be found in [321 ■ 
We have checked that, through the above \ 2 function, we can accurately reproduce the published 
CHOOZ limits (analysis A of 32 ) in the parameter space p = {Am 2 ,sf 3 }. For the sake of precision, 
in this work we have used the most general 3u parameter space for CHOOZ 85J 

p = {±AmVm 2 ,44} • (69) 

The effect of the subdominant parameters [5m 2 , s\ 2 , sign(±Am 2 )] is, however, rather small in the current 
data analysis. 



7.2 KamLAND 

The published KamLAND energy spectrum [2011*2] contains 258 events (background+signal), which 
we analyze through a maximum-likelihood approach [20], described in detail in [03] • In particular, the 
KamLAND (KL) x 2 function is defined as 

Xkl(p) = -2 In max £ K l(p, a, a', a") , (70) 

where a parametrizes a systematic energy-scale offset, while a' and a" represent free normalization 
factors for two (poorly constrained) background components The above likelihood function is 

factorized as 

£ K l(P, «. a 'i a ") = £rate(P, «, «") X £ sh ape(P, «, a") X £ syst (a) , (71) 

where the first, second, and third term embed the probability distribution for the total rate, for the 
spectrum shape, and for the systematic offset a, respectively; explicit expressions are reported in [93J. 
In particular, the spectrum shape term is further factorized into the probability distributions D(Ei) for 
finding the 258 KamLAND events with observed energies {-E'i}i=i,...,258 : 

258 

£shapc = II D (^) • ( 72 ) 
1=1 

A final remark is in order. In [03] the KamLAND analysis was performed in the 2v parameter space 
p = {5m 2 , Si 2 }, where the published bounds [201 121] have been accurately reproduced. In this work we 
have instead used the full 3z/ parameter space relevant for KamLAND, 

p = {5m 2 , s 2 12 , s 2 13 } . (73) 

We have checked, for a number of representative points in the 2>v parameter space, that the addition 
of KamLAND time-variation information [03] does not alter in any appreciable way the KamLAND 
constraints in such 3z/ space. 
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7.3 SK atmospheric data 



Our SK data analysis includes the zenith angle distributions of leptons induced by atmospheric neutri- 
nos, for a total of 55 bins, as discussed in |117j . With respect to |117j . we have updated the experimental 
event rates and their statistical errors { i?° xpt ±<r^ tat } from We also use three-dimensional neutrino 
fluxes j! 26| for the calculation of the theoretical rates i?J; heo (p). For convenience, we normalize both 
the experimental and theoretical rates to their no-oscillation value in each bin, as shown in Fig. 18. We 
consider eleven sources of systematic errors, which can produce a shift of the theoretical predictions 
through a set of "pulls" {£fc}fc=i,...,n |H7| . 

C e °(p)=< he °(p) + E^, (74) 

k=l 

where the response cjj of the n-th bin to the fc-th systematic source is numerically given in [117] . The x 2 
function is then obtained by minimization over the ^'s (which are partly correlated through a matrix 



xIk(p) = mi ? 

t?fei 



n=l \ u n I k,h=l 



(75) 



Minimization leads to a set of linear equations in the ^'s, which are solved numerically. The solution 
£ k can provide useful statistical information about the preferred systematic offsets and theoretical rate 
shifts, as discussed in Sec. 4.1. 

Finally, the 3z/ parameter space used for the SK data analysis in this work is 

p = {±Am 2 , S23, s 2 3 ,cos<5 = ±l,5m 2 *, s^} , (76) 

where 5m 2 * and si*, are the (fixed) best-fit values from the solar+KamLAND data analysis in Sec. 3.1. 
We have verified, in a number of representative points, that variations of these two parameters within 
their ±2cr limits do not alter the SK atmospheric data analysis in any appreciable way. 



7.4 K2K 

In this work, the K2K analysis is based on a 6-bin energy spectrum as in |117j . but including updated 
data [HH I129j as shown in Fig. 17. We cannot perform a K2K spectral analysis with finer binning (as 
the official one pH] ) for lack of published information about bin-to-bin systematic errors and their cor- 
relations in the last data release. The x 2 definition is based on a pull approach (with 7 systematic error 
sources |117j ). but the small number N of events in each bin requires a Poisson statistics, implemented 
through [TTII 

»theo v 



xLk(p) = mi V 



2 £ I nT° - ivr pt - ivr pt in ^ 1 + e ^fewa 

n=l \ iVn / k,h=l 



(77) 



with shifted predictions 

AC C °(p) = Ae eo (p) + E^. (78) 

k=l 

Numerical values for the K2K response functions and for the correlation matrix p s ^ are given in 

cm- 

For the sake of precision, and for consistency in the SK+K2K(+CHOOZ) combination, we have 
used for K2K the same 3z/ parameter as for SK 

p = {±Am 2 ,S23,s 2 3 ,cos5 = ±1, 5m 2 * , sf* 2 } , (79) 

although the subleading effects of the last two parameters are rather small in the K2K data analysis. 
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7.5 Solar neutrinos 



The definition of the solar neutrino x 2 is rather complex, both because it includes 119 observables and 
55 sources of systematics, and because it currently involves also correlations of statistical errors in the 
SNO salt data [17] and of systematic error sources in the BS05 (OP) solar model "9*8*1 19*9*1 H-UOj . Here 
we will mainly highlight the differences of the new x 2 inputs, with respect to our previous discussion in 

Ug. 12 

The formal x 2 definition is based on a pull approach, 



Xsoi(p) = mhi 

t?fe> 



119 55 

E X n[P~ ]nmX m + E £h[p7y B t]hk£k 
n,m h,k 



where 

*„ = <-(p)-^ , (81) 

and 

55 
k=l 

The presence of statistical error correlations (p nm ^ S nm ) in the latest SNO data [T7j does not spoil 
[93] the advantages of the pull approach discussed in jMj (where p nm = S nm ). We remind that the 
parameter space used in the solar neutrino data analysis is 

p = {<W,s 2 2 ,s 2 3 } , (83) 

except for Fig. 13 and related comments, where the very weak sensitivity to iAm 2 has been discussed. 

In the following we briefly describe, in ascending order, the n = 1...119 observables and the 
k — 1 . . . 55 systematic error sources of the solar neutrino analysis. 



7.5.1 Observables 

n — 1 ( Chlorine rate). The Chlorine rate input [j"J[""""] ^ s RT Pt ^ za i = 2.56±0.23 SNU. A p-dependent 
cross-section error is added in quadrature to <j[ (as described in [51]) to obtain the total uncorrelated 
error <j\. 

n — 2, 3 ( Gallium total rate and winter-summer asymmetry) . The Gallium (GALLEX/GNO+SAGE) 
input for the total rate is i?2 Xf>t = 68.1 ±3.75 SNU [01 HU HE]. A p-dependent cross-section er- 
ror is added in quadrature to o~' 2 (as described in "0""]) to obtain the total uncorrelated error The 
combined (GALLEX/GNO+SAGE) input for the winter-summer asymmetry (corrected for geometrical 
eccentricity effects [JTH]) is R 3 ±a 3 = -0.6 ± 7 SNU [TJ]. 

n = 4, ... ,47 (SK solar neutrino spectrum in energy and nadir angle). The analysis of the SK 
44-bin spectrum jTSj uses the same experimental input as described in appendix C of [HI], but all 
the theoretically computed quantities have been updated (in each point of the parameter space p) to 
account for the new hep and 8 B solar neutrino input in the BS05 (OP) model [0E1I1OO . 

n = 48, ... ,81 (SNO spectrum, pure D 2 phase). In this data set 45., events from NC, CC and 
ES scattering (and from various backgrounds) are not separated, and the global NC+CC+ES energy 
spectrum is analyzed. The spectrum information includes 34 bins (17 day + 17 night [33]) and is treated 
as described in appendix D of jHlj. As for the above (SK) data set, the only update concerns theoretical 
calculations, in order to account for BS05 (OP) solar model hep and 8 B input. 

12 Such new input includes, for all solar neutrino observables, also the electron density and neutrino production profiles 
from the BS05 (OP) SSM which are relevant for calculating solar matter effects on neutrino flavor evolution. 
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n = 82, . . . , 119 (SNO CC spectrum and ES+NC rates, salt phase). In this recent SNO data set |17j . 
the addition of salt has allowed a statistical separation of CC, NC and ES events. The CC spectrum 
includes 34 bins (17 day and 17 night, n = 82, . . . 115). Four more data points concern the total NC rate 
(day and night, n = 116, 117) and ES rate (day and night, n = 118, 119). The experimental values for 
the rates and their statistical errors (i?° x pt ± a n ) are taken from Table XXX (CC) and XXIV (NC+ES) 
of [T7j. These are the only data where the (statistical) errors are correlated, namely, p n ,m>&2 ^ 8 nm . 
The two (19 x 19) statistical correlation matrices for day and night data are taken from Tables XXXII 
and XXXIII of [Ej, respectively. 



7.5.2 Systematics 

k — 1, . . . , 20 (Standard Solar Model systematics). In jM] 12 sources X k of systematic errors were 
considered in the input solar model, individually denoted as Sn, S 33 , S34, S^u, Si 7, Lum, Opa, Diff, 
C~Be, Shep, and Z/X, each source being affected by a relative error AX k /X k (see also |16(Jj and refs. 
therein). With respect to [9*1] . the first 11 systematic uncertainties are unchanged, except for the 
updated value of AX/X for S 17 , currently set to 0.038 [M] (it was 0.106 in [Ml fijjft). The former 
systematic error source Z/X (solar metallicity) is now separated [HE] into 9 elemental uncertainties 
related to C, N, O, Ne, Mg, Si, S, Ar, and Fe, whose AX k /X k values are taken from the conservative 
estimate in Table 4 of [HE]- Three among these new error sources (O, Ne, Ar) are fully correlated [98J, 
and make the matrix p s ^ k st different from unity for (and only for) h,k — 14, 15, 19. The 11+9=20 SSM 
systematic error sources X k affect the neutrino fluxes through log-derivatives a ik (see 1160". |98j). 
For k = 1, . . . , 11, such derivatives are unchanged with respect to [HH I160j . For the new metallicity 
uncertainties (k = 12, ...,20), we take the log- derivatives from Table 1 of 98J. Finally, all SSM 
systematic uncertainties are propagated to the final c k n values in each point of the parameter space, as 
described in [9*"""| . 

k = 21 (^B spectrum shape uncertainty). The treatment of this systematic error |119j (which affect 
all solar neutrino observables) is unchanged from [9*4*j . 

k — 22, . . . , 32 (SK spectrum) and k — 33, . . . , 39 (SNO no-salt spectrum). The SK spectrum data 
{i? n } n=4 47 are affected by 11 systematics k = 22, ...,32. Analogously, the SNO (pure D 2 phase, 
{Rn\n=4s, ...,8i) data are affected by 7 systematics. They are treated as described in f9*4*j . 

k = 40, . . . , 55 (SNO salt-phase systematics). The recent SNO data in the salt phase {-R n }n=82,...,ii9 
are affected (in equal way during day and night) by 16 sources of systematic errors [T7j. The response 
function c\ is numerically given in Table XXXIV of [T7], in terms of (generally asymmetrical) upper 
and lower values c k+ and c^~. Such values cannot be incorporated at face value in a \ 2 function, 
which postulates symmetrical uncertainties. As mentioned in Sec. 3, in the pull approach we solve 
this problem through the prescription suggested in [§"""], namely, by shifting the central values of the 
observables (through the error half-difference) and by attaching symmetrized (half-sum) systematic 
errors as 

55 k+ _ „k- 55 k+ _i_ „k- 

Rn = Rn+ E n 9 n + " Z n (n= 82,...,119). (84) 

fc=40 1 fc=40 

7.6 Observables sensitive to absolute neutrino masses 

For the x 2 functions related to m^, m,pp, and S, we refer the reader to the thorough discussion given 
in [53]. With respect to [S3], in this work we have only reduced the theoretical nuclear uncertainty 
affecting mpp, according to the recent results reported in |154j (see Sec. 6.2). 
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Solar constraints at 2a (tf T3 = 0) 
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Figure 1: Regions separately allowed by the Chlorine (CI), Gallium (Ga), Super-Kamiokande 
(SK) and Sudbury Neutrino Observatory (SNO) experiments at the 2a level (A% 2 = 4) in the 
(5m 2 ,tan 2 # 12 ) plane, for 13 = 0. The LMA region allowed at 2a by the Cl+Ga+SK+SNO 
combination is superposed in each panel. 
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v e potential and production regions (BS05) 
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Figure 2: Neutrino potential V = y/2GpN e as a function of the normalized radius in the 
Sun. Also shown are the radial production regions for 8 B, 7 Be, and pp solar neutrinos (in 
arbitrary vertical scale). The curves refer to the Bahcall-Serenelli 2005 standard solar model. 
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Figure 3: The energy profile of the solar v e survival probability P ee for best-fit LMA values 
and #i3 = 0. The function P ee (E) shows a smooth transition from vacuum to matter- 
dominated regime as E increases, with some differences induced by averaging over different 
production regions (for 8 B, 7 Be and pp neutrinos) and, to a smaller extent, by nighttime (N) 
Earth effects with respect to daytime (D). Also shown are the corresponding solar neutrino 
energy spectra (in arbitrary vertical scale). 
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Figure 4: KamLAND constraints in the mass-mixing plane (5m 2 , sin 2 6 I 12 ) and for 13 = 0, 
as derived by an unbinned maximum-likelihood analysis of the total rate, spectrum shape, 
and rate+shape information. Contours are shown at 1, 2, and 3a level. 
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Figure 5: Solar and KamLAND constraints in the mass- mixing plane (5m 2 , sin 2 12 ) and for 
#13 = 0, shown both separately and in combination, at 1, 2, and 3a level. 
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Figure 6: Isolines of the charged-to-neutral current flux ratio and of the day- night asymmetry 
in SNO, superposed to the global (solar+KamLAND) LMA allowed region. 
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SK and SNO response functions 
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Figure 7: SNO response function for electron kinetic energy threshold T e > 5.5 MeV, and 
equalized SK response function at T e > 7.21 MeV. For such threshold, the SNO (CC) and 
SK (ES) solar v measurements are sensitive to the same parent neutrino energy spectrum, 
up to the negligible differences between the two curves in the figure. 
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Figure 8: Model-independent SK+SNO constraints on the flux of solar 8 B neutrinos and 
on their survival probability P ee (at 2a), averaged over the common response function in 
Fig. 7. Also shown is the ±2<r range of from the BS05 (OP) standard solar model. 
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Bounds on a MSW (Solar + CHOOZ + KamLAND) 
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Figure 9: Bounds on the amplitude of matter effects, taken as a free parameter czmsw- The 
cases omsw = 1 an d «msw — correspond to standard and no matter effects, respectively. 
The vertical axis represents the number of standard deviations (i.e., y/Ax 1 ) from the best- 
fit, as obtained by an analysis of all solar and reactor (KamLAND+CHOOZ) data with 
marginalized (5m 2 , sin 2 62) parameters. Standard matter effects are confirmed within a factor 
of ~ 2 uncertainty at 2a. 
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Figure 10: Pulls (at the global best-fit point for #13 = 0) of all the observables included 
the solar neutrino data analysis. 
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Figure 11: Pulls (at the global best-fit point) of the systematics included in the solar neutrino 
data analysis. 
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Figure 12: Preferred shifts of the solar neutrino fluxes (at the global best-fit point) with 
respect to the central values of the BS 2005 (OP) standard solar model. 
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neutrino survival probability from P 2u (#13 = 0) to P 3u (#13 > 0). See the text for details. 
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Figure 14: Three flavor analysis of solar and KamLAND data (both separately and in com- 
bination) in the parameter space (5m 2 , sin 2 612, sin 2 #13). The contours represent projections 
of the region allowed at 2a (A% 2 = 4). 
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Figure 15: Analysis of SK atmospheric neutrino data in the plane (Am 2 , sin 2 #23) at 13 = 0, 
for increasingly accurate inputs. The curves represent bounds at 1, 2 and 3a level. Upper 
panel: one-dimensional input fluxes and 5m 2 = 0. Middle panel: three-dimensional input 
fluxes and 5m 2 = 0. Lower panel: three-dimensional input fluxes and (5m 2 , sin 2 9u) fixed at 
their best-fit LMA values. 
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Figure 16: Analysis of SK and K2K data (both separately and in combination) in the plane 
(Am 2 , sin 2 ^23) at # 13 = 0. The parameters (5m 2 , sin 2 # 12 ) have been fixed at their best-fit 
LMA values. 
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Figure 17: K2K event spectrum as a function of reconstructed neutrino energy. Data (as used 
in this work) are shown by dots with ±ler statistical error bars. The histograms represent 
our calculations for no oscillation (dotted), and for oscillations at the SK+K2K best fit in 
Fig. 16 (dashed: with no systematic shifts; solid: with systematic shifts allowed). 
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Figure 18: SK distributions of leptonic events as a function of the cosine of the lepton zenith 
angle 9z, normalized to no-oscillation expectations in each bin. From left to right, the data 
samples refer to sub-GeV electrons (SGe), multi-GeV electrons (MGe), sub-GeV muons 
(SG/i), multi-GeV muons (MG/i), upward stopping muons (UP//), upward through-going 
muons (UT/i). Data are shown by dots with ±ler statistical error bars. The histograms 
represent our calculations at the SK+K2K best fit in Fig. 16 (dashed: with no systematic 
shifts; solid: with systematic shifts allowed). 
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Figure 19: Pull analysis (bin-by-bin) of the SK observables in Fig. 18. 
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Figure 20: Pull analysis of the eleven systematic error sources included in the SK analysis. 
From top to bottom, systematic pulls refer to: overall normalization (a); energy spectrum 
slope (5); /i/e flavor ratio for sub-GeV (f3 s ) and multi-GeV (/3 m ) events; relative normalization 
of partially and fully contained events (p); normalization of all upgoing muon events (p s ) 
and of through-going muon events alone {fit)] up-down asymmetry error of sub-GeV {rj s ) 
and multi-GeV {r] m ) events; horizontal/vertical error of low-energy {Kf) and high-energy 
(ac m ) events. 
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Figure 21: Representative examples of subleading three-neutrino effects in the SGe and MGe 
samples. See the text for details. 
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Figure 22: Isolines of the up-down electron asymmetry in the sub-GeV atmospheric neutrino 
sample, normalized to no-oscillation expectations. Left and right panels refer to normal and 
inverted hierarchy, respectively. Upper panels: 5m 2 = 0; middle panel: (5m 2 , sin 2 9 12 ) fixed 
at their best-fit values and 5 = 0; lower panel: as for the middle panel but with 5 — ir. 
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Figure 24: Three-neutrino analysis of SK+K2K+CHOOZ data, including subleading LMA 
effects. The results, marginalized with respect to (Am 2 ,^) and shown in terms of the 
X 2 (si3) function for both cos 5 = —1 (left half of the panel) and cos S = +1 (right part of the 
panel). The two CP-conserving cases cos 5 = ±1 smoothly merge at S13 = (vertical dotted 
line). The solid and dashed curves refer to normal and inverted hierarchy, respectively. 
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Figure 25: Three-neutrino analysis of SK+K2K+CHOOZ data, including subleading LMA 
effects. The results, marginalized with respect to s^ 3 , are shown in the (Am 2 , S23) plane for 
the two hierarchies (left and right) and the two CP-conserving cases (top and bottom). 
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Figure 26: Three-neutrino analysis of SK+K2K+CHOOZ data, including subleading LMA 
effects. The results are shown as projections of the (Am 2 , s% 3 , s 2 3 ) allowed regions (at 1, 2, 
and 3cr), marginalized with respect to the four cases [cos5 = ±1] <S> [sign(±Am 2 ) = ±1]. 



58 




sin 

Figure 27: Global three-neutrino analysis of oscillation data. Bounds on s 2 3 are shown 
in terms of na = ^Ax 2 for KamLAND (dotted curve), solar (dot-dashed curve), so- 
lar+KamLAND (short-dashed curve), SK+K2K+CH00Z (long-dashed curve) and all data 
combined (solid curve). In each case, the continuous parameters (Am 2 , s^-, s 2 3 ) and — if 
applicable — the discrete parameters [cos 5 = ±1] <8> [sign(±Am 2 ) = ±1] are marginalized 
away. 
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Figure 28: Global three- neutrino analysis of oscillation data. Bounds on each of the param- 
eters 5m 2 , Am 2 , sin 2 9i2, and sin 2 #23 are shown in terms of na = y/Ax 2 - In each plot, all 
parameters but the one in abscissa are marginalized away. 
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Figure 29: Regions allowed at 2a by the global analysis of neutrino oscillation data in the 
parameter space of non-oscillatory observables (m^, rripp, E). The regions are projected onto 
the three coordinate planes for both normal hierarchy (thick curves) and inverted hierarchy 
(thin curves). 
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Figure 30: Analysis of oscillatory and non-oscillatory observables in the plane {nipp, E). The 
2cr horizontal band is preferred by the positive 0^2/? claim, while the slanted 2a regions 
below are preferred by all other data (oscillation and cosmological data). The absence of 
overlap indicates tension among the data in the sub-eV range. 
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Figure 31: Analysis of oscillatory and non-oscillatory observables in the parameter space 
(171/3,171/313, E), including all data but the claimed lower bound from Oi/2/? searches. 
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Figure 32: Analysis of oscillatory and non-oscillatory observables in the parameter space 
(mp, 171/3/3, S), including the claimed Oz/2/3 signal but excluding the recent Lya forest data. 
The contours of the allowed regions for normal and inverted spectrum hierarchy overlap 
(regime of degenerate spectrum). 
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